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A correlation analysis of a thickness tendency equation including 
heating and vertical motion terns was carried out. One of the main pur- 
poses of the investigation is to provide a firmer foundation for a 1000- 
500 mb thickness forecast scheme. 

A large data sample was used and single and multiple linear coeffi- 
cients were calculated. These were computed for a wintertime situation 
with separate analyses for Atlantic and Pacific areas. 

The resulting correlation coefficients for the advection, vertical 
motion and heating terms gave substance to the hypothesis that the observed 
change in thickness may be fairly well approximated by the change due to 
horizontal advection only. The vertical motion term gave a smaller but 
significant correlation, while heat exchange terms did not correlate 
significantly. 
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1 . 



INTRODUCTION 



The ability to produce analyses and prognoses of weather maps by 
numerical methods using high speed digital computers on an operational 
basis is now a reality. Most of the efforts in numerical weather pre- 
diction have been in relation to the pressure and temperature fields; 
however some pioneering attempts have been made to predict precipitation, 
clouds and visibility. 

Perhaps the most successful forecasts have been of the 500-mb 
contours. However, somewhat less success has been achieved in obtaining 
a good surface or 1000-mb prognosis. Many schemes have been tried with 
varying degrees of success. One method that has a great deal of promise 
makes use of the 500-mb and 1000-500 mb thickness prognoses to obtain a 
1000-mb forecast. The U. S. Naval Fleet Numerical Weather Facility 
(FNWF) has utilized an empirically adjusted advective technique for 
thickness prediction. A more complete thickness forecast model includ- 
ing vertical motion, heating and friction is presently under investiga- 
tion by Professor G. J. Haltiner of the U. S. Naval Postgraduate School. 
This program utilizes a thickness tendency derived from the first law of 
thermodynamics and is shown in a subsequent section of this paper. In 
some of the earlier applications of this equation some of the terms have 
been omitted by way of simplification because of the lack of information 
necessary to calculate these terms, for example, diabatic heating. Re- 
cently FNWF has made computations of the heat exchange on a regular 
basis. With this data now available it is appropriate to consider 
whether a better estimate of a thickness tendency might be obtained if 
more terms of the equation were retained and fewer simplifying assumptions 
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made . 

The purpose of this investigation is to conduct a first step in a 
determination of the usefulness of this equation, that is, to statisti- 
cally correlate the observed 1000-500 mb thickness changes with estimates 
of the various terms in the thermodynamic equation. As shown by Panofsky 
and Brier [8] linear correlation coefficients give a measure of associa- 
tion between two variables which does not depend upon any arbitrary choice 
of units by which the original variables were measured. From the correla- 
tion coefficients obtained here each term will be examined as to its re- 
lation to the observed thickness change and its significance as a pre- 
dictor. From these results the value of this equation in a prediction 
scheme can be estimated. 
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2 . 



AREA OF INVESTIGATION 



FNWF routinely analyzes data over the Northern Hemisphere on a 
63 x 63 grid, with a grid distance of 381 KM. 

As shown later in this paper the scheme used in obtaining upwind 
points for data extraction was designed for use on a 22 x 22 grid extract- 
ed from any section of FNWF's 63 x 63 grid. The total heat exchange 
calculations were originally derived for use over the oceans only, how- 
ever FNWF calculates this term over the whole grid, and then plots only 
the ocean points. Because the land area points were to be kept at a 
minimum, all available values of the heat exchange were used. The known 
decrease in the accuracy of the value over land as compared to that over 
the ocean nevertheless was accepted in this study and recognized as a 
possible source of error. 

Due to the facts mentioned above two different areas of 22 x 22 grid 
size were selected for this study. These areas were chosen because most 
of the grid is oceanic with a minimum of continental area. 

Area 1, which will be referred to as the Atlantic area, is shown in 
figure 1. This area includes the North Pole and the Gulf Stream and has 
a minimum of equatorial area; the latter is near the boundary of. FNWF 's 
grid where results are considered to be less accurate. 

Area 2, which will be referred to as the Pacific area, is shown in 
figure 2. This area is almost entirely oceanic, but does contain more 
of the equatorial or boundary area. These areas will be treated separ- 
ately and the correlation coefficients will be calculated independently. 



11 



3. 



BACKGROUND 



The thickness prediction equation may be derived in the following 

way. The first law of thermodynamics can be written 

Q=C,(| tV'7T 1 

i ii iii iv v (1) 



Combining the equation of state and the hydrostatic equation gives 






(2) 



Substituting (2) into (1) in terms (ii) and (iii) gives 





O0 



(h _RT \ 

cff) 



(3) 



Rearranging, we have 

^) = . V . V (^) <*) 

Integrating between 1000 and 500 mb 

j^=-yVv(^U/> +fo-uUf, - f <Lp (5) 

where 

r-fr (-£-■£> 

Now integrating over a 12-hour time period gives 



Ah. = f) V- Wj^y/ 9 ' 6 * + ft <no^f 

0 /... _ <>U _ 



JLt 




(7) 
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The purpose of this investigation is to correlate the 12-hour thick- 
ness changes with numerical approximations of the terms constituting the 
right hand side of equation 7. Thus Ah. will be correlated with 
estimated mean values of the horizontal advection, A, the vertical motion 
term, W, and the total heat exchange, H, for the 12-hour period. This 
relationship may be written 

Ak = ^ W + H (8) 

The coefficients of the terms in equation 7 can be considered as 
merely constants of proportionality and thus can be omitted in the 
statistical evaluation. 

The parameter CT as defined here differs slightly from the value 
of CT” used by FNWF. However, here again the difference can be con- 
sidered essentially as a constant of proportionality and will not affect 
this study. 



13 



4. 



DATA COLLECTION 



The basic data for this study was obtained for the three-day period 
9 to 11 February 1966 from the Fleet Numerical Weather Facility. 

The initial fields utilized were the SD-500 field which is the dis- 
turbance component (perturbation) of the 500-mb height field, the SR-500 
field which is the residual component (basic flow) of the 500-mb height 
field, the SD-1000 field which is the disturbance component at the 1000 
mb level and the SR- 1000 field which is the basic flow at the 1000 -mb 
level . 

These fields were used in a routine originally designed for the 
interpretation of meteorological satellite observations as developed by 
Nagle, Clark and Roll [7] and modified by Clark for use in obtaining the 
advective change in thickness, A. The scheme entails generation of 
stream function fields, <s f // , figures 3 and 4, from the SD and SR fields 
averaged over the 1000-500 mb layer. By a timewise interpolation of the 
12-hour distributions, hourly distributions of ^ were obtained which 
would act as evolving frame- of- reference steering fields. From these 
V fields a horizontal velocity distribution was computed. This wind 
was based on the geostophic approximation and in this case was the mean 
wind in the layer. 

At this time the 1 2-hour ly winds for a 22 x 22 subsection of the 
63 x 63 grid were extracted. This subregion was then interpolated in 
both time and space to give hourly winds into an expanded 44 x 44 array 
based on the 22 x 22 section. This maneuver quadruples the number of 
data points contained in the region. 

By utilizing the stream function and the wind field, an upstream 
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trajectory was projected and the position of a column of air at time 
t - 12 was determined. This was done for each data point in the region. 

Checks in the program flagged any point whose upstream position was out 
of the 44 x 44 grid. This point was then discarded for correlation 
purposes. 

Several steps were involved in the procedure used to determine X. 
Thickness fields, figures 5 and 6, were obtained from FNWF's D-500 and 
D 1000 fields. These D fields give the difference between the height 
of the surface (500 and 1000 mb, respectively) and its corresponding 
standard height. The trajectory scheme was then applied to this field 
to yield a new parameter, 6, which is the combined local change of thick- 
ness, Alro , and the estimated horizontal advection, A. Figures 7 and 
8 depict this parameter in the Atlantic and Pacific areas, respectively. 
The strong gradient in the extreme western section seen in these and 
subsequent maps was due to the check-flag system discussed previously. 
These are the areas of inflow where the upstream position was off the grid 
causing a zero to be inserted. While the zero value was eliminated be- 
fore correlations were made, the point was plotted as such resulting in 
tight gradients where none actually exist. 

The local or observed 12-hour change in thickness, figures 9 and 
10, was produced by subtracting the field at time t - 12 from the field 
at time t Q . To get an estimate of A at each point, A ho was subtracted 
from h; this was not completed until the format of the data had been 
changed and the correlation program was in use. 

A further modification of the routine allowed W to be computed. 

This vertical motion term was obtained by taking the product of 0“ and 
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Uj , the vertical motion in (x,y,p,t) coordinates, and averaging in 
both time and space. 

To accomplish this sigma-1 and sigma-3, the static stability for the 
layers, 1000-750 and 600-750 mb respectively, weie utilized. These C” 1 s 
were derived from a model by Roll, Bibbo and Clark [2] which is charact- 
erized by pressure levels defining a number of layers within which the 
static stability is independent of pressure. 

These terms were multiplied by the c*J fields at 850 and 500 mb. 

FNWF computes these parameters by a numerical solution of a quasi - 
geostrophic omega equation, as shown by Haltiner, Clarke and Lawniczak 

[1]. The resulting fields were <TuJ ocri » figures 11 and 12, and , 

ojU jUU 

figures 13 and 14. The trajectory scheme was again used to obtain up- 
wind values of these terms at time t - 12 which were then averaged with 

o 

values at t Q to yield the 12-hour averages, cf^O and shown 

in figures 15, 16, 17 and 18. The final averaging of and 

to obtain W was completed by the correlation program itself. 

The same procedure was used to obtain H. This term was derived from 
the fields of Q^, the total air-ocean interchange as produced by FNWF 
shown in figures 19 and 20. The computation program of the total heat 
exchange by numerical methods and its adaption to the computer evolved 
from earlier manual methods developed by Laevastu [4]. The heating func- 
tion was obtained from the following equation: 

- Qs ~ Q.r ~ Qb “ Qfv. “ Qe (9> 

-2 -1 

where all Q's will be expressed in units of g-cal cm (24 hr.) 

Because of the relative newness of the computerized form of this product 
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each term will be briefly discussed. 



The insolation, Qs, is calculated from the equation 



Q s = o. ih- A k ^ ( \ ~ O. OOO&C ) 



( 10 ) 



which takes into account the total time that the sun is above the horizon, 

t., the noon altitude of the sun, A , and a cloudiness factor, C. Since 
d n 

it is virtually impossible to take into consideration various types of 
clouds and their different effects on the insolation, Laevastu [3] 
devised an average cloudiness factor which requires only a knowledge of 
total cloud tenths. The results give a reasonable estimate of the total 
incoming radiation. 

The albedo, Q^, is an estimate of the radiation reflected from the 
surface. Laevastu [5] statistically determined an empirical formula for 
calculating this term. 



This formula is valid only for daily calculations. A simplified formula 
for the computation of short term (hourly or 3-hourly) albedo is 



The effective back radiation, Q^, is the long wave radiation given 
off by the sea surface. This represents a heat loss at the surface 



and makes use of a linear formula of Lonnquist [6] which is corrected by 
Moellers cloudiness factor. 




(ID 



Qc “ Q s = Average solar altitude 



( 12 ) 



Q b = (zw - I.S6 Tv - 0.9fUo)( I - o-ous-c) (13) 
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The convective transfer of sensible heat, Q, , was obtained from 

h 

Qw = 2 ^ ( O. + 0.11 \/)(T w - 7 Sl) j 7^ > 71 (14) 

when heat is transferred from the water to the air; and 

* 3V(T» - To.) ( T w < T«. (15) 

when heat is transferred from the air to the water. 

The latent heat transfer, Q , is that heat addition or loss due to 

e 

evaporation or condensation. Lavastu [5] shows that a modified formula 
of Rohwer has been found to be most accurate in estimating evaporation; 
namely, 

Qe. * (o. 2 t + o .->7 v)(o.?8 if o.qae w >e«, <i6) 

A modified form of this equation can be used for condensation at the sea 
surface, 

Q c - 0.O77 v(o.iee„ - e A ) L t f rf o,9f^<e^ ( i 7 ) 

To calculate these terms and get a total heat exchange product 
certain basic synoptic numerical analyses must be available. These are: 
Surface winds, surface air temperature and water vapor pressure and 
clouds. These are all routinely produced at FNWF. 

By utilizing the trajectory program again, the t Q - 12 values of 
Q were obtained. Averaging over the 12 hours produced H, figures 21 
and 22. No further data was needed. 
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5. 



OBJECT AND METHOD INVESTIGATION 



The purpose of this investigation was two-fold. An attempt was 
made a) to statistically correlate the observed change in thickness 
with estimated values of the terms on the right hand side of equation 
(7) and b) to estimate the significance of each term in a prediction 
scheme . 

To accomplish this a multiple regression and correlation analysis, 
the BMD03R programmed for the CDC 1604 electronic computer and converted 
for COOP MONITER Fortran 63, was utilized. This program will perform 
linear-multiple correlation analysis on data in any combination by allow- 
ing the designation of a dependent variable and deletion of sets of in- 
dependent variables. This can be done as many times as desired. 

At the completion of the data collection phase there were five 
parameters stored for each data point over a total of 8229 Atlantic and 
6739 Pacific points which were considered good for correlation purposes. 
These variables were: (1) tx ; (2)oraJ g50 ; (3)cru' 500 ; ^ ** » and (5)Ak. 
To place this data in a form appropriate for the BMD program the format 
had to be changed from a row-wise listing to a column presentation. 

After accomplishing this, spot checks were made to ensure that the data 
had not gotten out of phase during the transformation. 

Since the parameters A and W had not actually been obtained as yet, 
the transgeneration feature of BMD03R was used. This allows for the 
generation of new variables by performing certain mathematical opera- 
tions on the original variables. Here A, variable 6, was obtained by 
the subtraction (1) - (5) = (6) and W, variable 7, was produced from 
[(2) + (3) ]% = (7). The data was then ready for the correlation analy- 
sis . 
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6 . 



RESULTS AND DISCUSSION 



The data was divided into Atlantic and Pacific regions with identi- 
cal calculations made for each area. Initially, correlation analyses 
were carried out to show relationships between pairs of variates. Many 
combinations of independent and dependent variables were used. The re- 
sults are shown in Tables 1 and 2. 

After examining the figures, multiple linear correlations were made. 
This was done by progressively adding the term from equation (8) with the 
highest partial correlation coefficient. 

As a side operation, a multiple correlation analysis was made on 
the terms of the equation 

k oC W “V* [\ 

These results are shown in Tables 3 and 4. 

The values of R obtained from the Atlantic data were slightly higher 
than those in the Pacific in almost every case. This can be partially 
explained by a) the scarcity of reporting stations in the Pacific result- 
ing in larger interpolated areas, and b) the Pacific area contains a 
larger band of equatorial or boundary region where the analysis is known 
to be less accurate. 

The coefficient obtained in the linear correlation of and H was 
extremely low in both the Atlantic and Pacific. 

The coefficient resulting from the correlation of /ikfvand W was lower 
than expected but was approximately the same in both Atlantic and Pacific. 

^ was correlated with 0~kJ oirr . and <TCJ separately to see if it was 

850 500 

necessary to obtain W or if the value of the parameter obtained at a 
single level would be sufficient. In the Atlantic there is little 
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Dependent 


Independent 


R 


R 2 


Variable 


Variable 






jCS la- 




.8599 


.7394 


Ah. 


w 


.1936 


.0375 


Ah- 




.0655 


.0043 


Ah. 




.1729 


.0299 


A h. 


O-i^s-oo 


.1850 


.0342 


W 


"A 


.4513 


.2037 


v7 


TT 


.2200 


.0484 


“A 


TT 


.1799 


.0324 




TT 


.2425 


.0588 


Irv. 


W 


.5802 


.3366 



Table 1. Results of non-multiple linear correlation analyses, 
Atlantic. R is the correlation coefficient and R^ is the co- 
efficient of determination. 
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Dependent 


Independent 


R 


R 2 


Variable 


Variable 






Ah- 


A 


.7967 


.6347 


Ah- 


W 


.1752 


.0302 


AK- 


PT 


.0985 


.0097 


AVv 




.2440 


.0595 


Ah, 


CHJsao 


.0825 


.0068 


v7 


A 


.2569 


.0660 


W 


“H 


.0404 


.0016 


A" 


TT 


.0962 


.0093 


K, 


TT 


.0483 


.0023 


FL 


W 


.5373 


.2886 


Table 2. 


Results of non-multiple 


linear correlation 


analyses, 



Pacific . 
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Dependent 


Independent 


R 


R Z 


Variable 


Variable 






AW. 


"A 


.8599 


.7394 


AW, 


"A + W 


.8871 


.7869 


AW. 


"A +W+TT 


.8889 


.7901 


H 


vs/ + ~VT 


.5920 


.3505 


Table 3. 


Results of the multiple-linear 
Atlantic. 


correlation 


analysis. 


Dependent 


Independent 


R 


R 2 


Variable 


Variable 






AK. 


"a" 


.7967 


.6347 


A W. 


A +W 


.8562 


.7331 


AW. 




.8625 


.7439 


H 


W + TT 


.5418 


.2936 


Table 4. 


Results of the multiple-linear 
Paci f ic 


correlation 


analysis. 
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difference in the coefficients but in the Pacific, the three differ 
considerably. 

The multiple correlation coefficient in the Atlantic increased 
from .8599 to .8871 to .8889 as W and H were added to "A. The increase 
in the Pacific was somewhat greater, going from .7967 to .8562 to .8625 
as these terms were added. Here the inclusion of the vertical motion is 
more significant. 

The analysis of with W and H was attempted with the thought 
that K. might be used in a cloud forecast scheme currently under study. 
Although the correlations in both areas are significant, they are certain 
ly not high enough to consider replacing vertical velocities by K. in a 
cloud forecasting scheme. 

The remaining analyses clearly show no significance and will not be 
discussed further. 
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7. 



CONCLUSIONS AND RECOMMENDATIONS 



As a result of this study, which was made using wintertime data 
only, the following conclusions were reached. 

(1) Observed changes in thickness over a 12-hour period are due 
almost entirely to horizontal advection, at least as approximated in this 
study. This result lends strong support to the simple advective scheme 
for predicting thickness which has been used operationally by FNWF dur- 
ing the past several years, i.e. 

Ah. <c V* V h. 

(2) The vertical motion term, W , though of secondary importance, 
does contribute significantly to the change in thickness and should be 
retained in a prediction scheme for increased accuracy. Of course, a 
decision must be made as to whether the improvement in forecasting skill 
is worth the time and effort necessary to include vertical motion. 

(3) The heat exchange term, H, as presently being calculated at 
FNWF and as utilized in this study contributes virtually nothing to the 
change in thickness and should be omitted in a thickness prediction scheme. 

As pointed out earlier, this study dealt only with winter data. An 
obvious and necessary extension should be the application to data from 
the other seasons, especially summer. 

It is hoped that the results of this study will provide guidance to 
individuals or groups engaged in the development of thickness prediction 
schemes. 
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Figure 1 Area of Investigation, Atlantic 




Figure 2 Area of Investigation, Pacific 
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Figure 5 1000 - 500 mb thickness field, h, Atlantic 
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Figure 6 Thickness field, h. Pacific 
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Figure 10 Local change of thickness, Ah, Pacific 




Figure 11 CTU) Atlantic 




Figure 12 (T6J 85Q , Pacific 







Figure 13 cn5c, 5oo* Atlantic 
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Figure 14 C TOJ ^qq* Pacific 
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Figure 16 ^^gSO’ Pacific 




Figure 19 Total heat exchange, Q , Atlantic 




Figure 20 ' Total heat exchange, Q , Pacific 
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Figure 22 Mean. Heating term, H, Pacific 
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